Table 1. Line equivalent widths (A) and veilings. 



Star UT Date Ha H/3 H7 HS He I Na Fe I Li I Ca II Ca II Ca K Fe I Fo II Fe II Veiling 



red/blue D1+D2 A6192 A6708 A8662 A8498 A5497 A4924 A4352 red blue 



AA Tau 


1989 Oct 13/13 


14.2 


2.2: 


3.5: 


2.7 


1.2 


-4: 


-0.3 


-0 


6 


0.5 


0.7 


2.2 


-0.3 


X 


X 


0.3 


0.2 


AS 205 


1987 Jun 10/- ■ ■ 


111.1 








1.6 


2.5 


-0.2 


-0 


1 


20.4 


39.7 




-0.1 






0.2 




AS 209 


1987 Oct 06/06 


113.4 33.5 


11.5 


7.9 


3.3 


2.9 


-0.3 


-0 


4 


16.6 


21.0 


21.6 


-0.2 


X 


X 


0.3 


0.3 


AS 353 


1987 Sep 06/06 


56.0 


12.6 


5.9: 


-2.9 


2.6 


2.9 


0.7 


-0 


1 


37.9 


36: 


X 


0.0 


5.6 


3.8 


1.7 


2.4 


BP Tau 


1986 Dec 21/21 


48.6 


10.7 


4.5 


4.8 


1.0 


-1.1 


-0.2 


-0 


4 


4.1 


3.5 


11.6 


-0.2 


X 


X 


0.5 


0.9 


CI Tau 


1987 Oct 12/- ■ ■ 


110.1 








3.4 


-0.5 


-0.3 


-0 


4 


19.7 


19: 




-0.2 


3: 




0.2 




CO Ori 


1992 Nov 15 


2.9 


-2.2 


-3: 


-2.9 


-0.1: 


-2.2 


-0.2 


-0 


2 


X 


0.2 


X 


-0.2 


X 


X 


0.2 


0.2 


CW Tau 


1998 Nov 03 


109.612.3 


X 


X 


1.9 


X 


-0.1 


-0 


1 


21.5 


28.3 


3: 


X 


2.2 


X 


1.5 


X 


DE Tau 


1988 Feb 04/04 


66.5 


44.0 


34.6 


X 


2.6: 


X 


X 


-0 


4 


7.1 


7.9 


76.4 


X 


X 


1.7 


X 


X 


DF Tau 


1992 Nov 14 


61.0 


6.8 


3.0 


3.2 


3.0 


0.8 


X 


-0 


2 


4.9 


5.3 


7.5 


X 


0.8: 


0.1: 


1.8 


1.8 


DC Tau 


1993 Oct 30 


60.5 


17.7 


8.2 


5.4 


2.4 


2.5 


0.3 


-0 


2 


41.4 


44.2 


20.8 


-0.1 


5.0 


3.0 


2: 


2: 


DK Tau 


1987 Oct 12/12 


16.7 


4.4 


0.9 


2.1 


1.0 


-2.1 


-0.2 


-0 


4 


2.1 


2.7 


3.5 


-0.2 


X 


X 


0.4 


0.7 


DL Tau 


1998 Nov 03 


110.2 24.8 


7.2 


3.4 


5.4 


7.0 


-0.1 


-0 


2 


39.9 


46.5 


8.5 


X 


7.3 


X 


1.9 


X 


DN Tau 


1986 Nov 12/12 


15.7 


5.6 


1.9 


3.9 


0.3 


-3.4 


-0.3 


-0 


6 


0.4 


0.5 


10.6 


-0.3 


X 


X 


0.1 


0.2 


DO Tau 


1998 Nov 03 


58.2 


5.3 


X 


X 


2.1 


2.4 


X 


-0 


2 


8.8 


7.9 


0.9: 


X 


1.9 


X 


0.9 


X 


DQ Tau 


1993 Dec 24 


157.123.4 


3.3 


2: 


2.7 


1.8 


-0.1 


-0 


3 


2.7 


2.3 


X 


X 


2.0 


X 


1.5 


1.5 


DR Tail 


1 QQ8 Nov 03 


81.0 


5.3 


1.0 


1.5 


4.2 


4.0 




-0 


I 


23.3 


25.4 


11.2 




2.4 


0.7 


2.5 




DS Tau 


1987 Oct 11/11 


27.3 


7.2 


4.6 


3.7 


0.9 


-0.5: 


-0.2 


-0 


2 


0.4 


0.3 


19.9 


-0.1 


X 


X 


0.6 


0.5 


GO Tau 


1989 Jan 17/18 


40.2 


8: 


4: 


2: 


0.8 


-2: 


-0.2 


-0 


5 


2.3 


2.4 


X 


-0.2 


X 


X 


0.3 


1.1 


GK Tau 


1988 Feb 04/04 


34.0 


9.6 


6: 




1: 


X 


X 


-0 


6 


4: 


1.1 


6.9 


-0.3 


X 


X 


X 


0.1 


CM Aur 


1987 Oct 11/11 


63.5 


12.4 


4: 


4.8 


0.7 


-2.1 


-0.3 


-0 


4 


0.3 


0.4 


32.5 


-0.3 


X 


X 


0.2 


0.1 


GW Ori 


1993 Dec 23 


31.5 


2.9 


X 


X 


X 


1.4 


-0.2 


-0 


2 


4.7 


5.2 


X 


-0.1 


X 


X 


0.2 


X 


HL Tau 


1986 Nov 12/- ■ ■ 


34.7 








1.2 


3.8 


-0.2: 


-0 


2 


28.2 


30.6 




X 


2.5: 




0.9 




RW Aur 


1986 Dec 21/21 


45.9 


4.1 


-2.1 


-1.2 


0.6: 


-8.3 


0.6 


-0 


2 


27.5 


37: 


14.0 


X 


X 


X 


2.0 


2.9 


RY Tau 


1986 Dec 21/21 


12.4 


-1.4: 


-0.8: 


-0.9 


-0.3 


0.2 


-0.2: 


-0 


2 


4.7 


5.2 


2.2 


-0.2 


X 


X 


0.1 


0.0 


SU Aur 


1986 Nov 14/14 


6.1 


-1.7 


-1.1: 


-1.6: 


-0.3 


-1.0 


-0.2 


-0 


2 


-1.0 


-0.7: 


X 


-0.3 


X 


X 


0.0 


0.0 


T Tau 


1989 Oct 13/13 


63.1 


9.8 


3: 


3.3 


0.5 


0.7: 


-0.2 


-0 


4 


11.2 


13.7 


27.7 


-0.2 


1.4: 


X 


0.1 


0.1 


UX Tau 


1986 Dec 22/- ■ ■ 


3.5 








X 


-2.0 


-0.3 


-0 


4 


0.2 


0.2 




-0.3 






0.0 




UY Aur 


1988 Feb 04/04 


97.0 


20.8 


10.7 




2.6 


X 


X 


-0 


4 


7.8 


12: 


96.2 


X 


X 


X 


0.1 


0.8 


V410 Tau 


1992 Nov 14 


2.9 


X 


X 


X 


X 


-3.4 


-0.3 


-0 


5 


0.2 


0.2 


2.2 


-0.3: 


X 


-1: 


0.1 


0.2 



Note. — When the lines could not be reliably used, the equivalent width and veiling were not measured (x) and when we 
did not have a spectrum at a given wavelength it is labeled as 
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Table 1. Velocity centroid (km s ^) of the line profile components. 





star 




Ha 






H/3 






H7 






H<5 




He I 


o : 




cl 


c2 


c3 


cl 


c2 


c3 


cl 


c2 


c3 


cl 


c2 


c3 


IN 




o ■ 

o ■ 
<N : 
































AA Tau 


-30.1* 


-34.1 


X 


-23.9* 


-26.9 


87.9* 


-1.4* 


-13.1 


103.5* 


-3.2* 


0.3 


101.0* 


-0.4 


9.0 


a ■ 


AS 205 


-109:* 


-46.4 


174.5 




















X 


-14.3 


■ 


AS 209 


-4.7 


X 


X 


1.3: 


X 


X 


8.1 


X 


X 


-1.4 


X 


X 


-4.2 


0.6 


^ ; 
































oo ■ 
T— 1 1 


AS 353 


-137:* 


-266:* 


30.1 


-179.2* 


97.7 


X 


-186:* 


131.6 


X 


-222.5* 


136: 


X 


X 


-41.2 


BP Tau 


4.4 


48.1* 


X 


-6.1 


63.8* 


X 


-0.9 


63.8* 


X 


-4.5 


63.4* 


X 


0.4 


-10.5 




CI Tau 


-98.3* 


-9.2 


X 




















-5.1 


-23.0 


> ■ 


CO Ori 


-37.0* 


-111.6* 


22.7 


-38.1* 


X 


X 


-38.4* 


115:* 


X 


-19.9* 


X 


X 


6:* 


X 


(N : 

(N . 
m ■ 


CW Tau 


-14.3* 


-130.9* 


-0.8 


-1.2* 


8.3 


X 


X 


X 


X 


X 


X 


X 


X 


-58.0 


DE Tau 


-0.4 


X 


X 


12.3 


X 


X 


8.2 


X 


X 








-2: 


-53: 


^ ■ 

o : 


DF Tau 


-70.8* 


-10.4* 


-5.3 


-70.2* 


-14.2* 


-17.9 


-65.3* 


-9.4* 


-28.3 


-57.3* 


1.0* 


-22.3 


4.6 


21.8 


o . 


DC Tau 


-102.9* 


-162.2* 


-38.4 


-150.1* 


-50.1 


X 


17.7 


X 


X 


21.0 


X 


X 


X 


-4.2 


o ■ 


DK Tau 


-26.3* 


-13.7 


X 


0.8* 


-14.3 


213.8* 


0.1* 


-7.8 


234.1* 


14.2* 


-30.6 


213.6* 


-0.2 


0.4 


ro-ph 


DL Tau 


-122.0* 


2.1 


59.1* 


-121.4* 


8.3 


X 


14.7 


X 


X 


21.0 


X 


X 


4.2 


-2.5 


DN Tau 


-58.2* 


-38.8 


X 


-43.1 


-7.7 


149.3* 


-2.4 


X 


X 


0.9 


-3.2 


X 


-2.8 


X 


DO Tau 


-76.0* 


-3.3 


43.1* 


-72.1* 


-6.3 


50.7* 


X 


X 


X 


X 


X 


X 


-0.5 


5.6 


































Civ: as 


DQ Tau 


-250.0* 


-8.2 


262.9 


-7.2 


297. 


X 


-6.7 


X 


X 


-12.1 


X 


X 


2.3 


-4.2 


DR Tau 


-155.0* 


-251.8* 


100.9 


-66.8* 


-194.0* 


100.6 


-44.7* 


-133.8* 


81.2 


-69.3* 


-125.8* 


58.3 


0.1 


-14.5 


DS Tau 


-54.8* 


-35.4 


X 


-62.1* 


-49.4 


X 


-65.7* 


-28.3 


X 


-33.9 


X 


X 


-7.4 


6.0 


% 


GG Tau 


-94.8* 


14.9 


57:* 


6.1 


X 


X 


10: 


X 


X 


42.2 


X 


X 


3.0 


-5.1 




GK Tau 


-28.3* 


-33.3 


X 


-24.4* 


-32.0 


X 


X 


X 


X 








-4.7 


X 




CM Aur 


-51.8* 


7.5 


X 


-43.7* 


-49.7 


X 


-45.9* 


-33.4 


249.3* 


-40: 


189:* 


X 


2.4 


-17: 




GW Ori 


-78.7* 


-5.2 


X 


-74.4* 


-22.6 


X 


X 


X 


X 


X 


X 


X 


X 


X 




HL Tau 


-38.7* 


-2.3 


X 




















X 


X 




RW Aur 


-96.0* 


4.5 


5.7* 


-88.7* 


9.5 


21.6* 


-76.8* 


34.9* 


X 


-84.6* 


-72: 


21.8* 


0.5 


-15.0 




RY Tau 


-12.2* 


14.3 


118.* 


28:* 


X 


X 


31:* 


X 


X 


36.7* 


X 


X 


54:* 


X 




SU Aur 


-22.4 


10.7 


103.1* 


38.9* 


X 


X 


20.4* 


X 


X 


23.6* 


X 


X 


29.2* 


X 




T Tau 


-171.* 


-2.68 


X 


-4.3 


X 


X 


9: 


X 


X 


2.7 


X 


X 


X 


-9.3 




UX Tau 


-48.8 


42.5* 


269.4* 




















X 


X 




UY Aur 


-30.6 


5.4 


X 


-15.0 


-21.9 


X 


-18.3 


-7.7 


X 








0.1 


10.1 




V410 Tau 


60.7 


102:* 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 


X 



Note. — Lines identified as "■ • ■" were not available and as "x" were not present or could not be reliably 
measured in the observed spectrum. Velocities marked with "*" represent absorption components, cl, c2 and 
c3 represent the emission and absorption components of the lines that could not be decomposed as NC and 
BO. 
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Table 1. Continued 





Star 


Na Dl 




Na D2 




Fe I 


Li I 




Call 




Call 




Ca K 


Fc I 


Fell 


Fell 


















A6192 


A6707 




A8662 




A8498 








A5497 


A4924 


A4352 


2000 




cl 


c2 


c3 


cl 


c2 


c3 


cl 


cl 


NC 


BC 


c3 


NC 


BC 


c3 


NC 


BC 


cl 


cl 


cl 


AA Tau 


4.3* 


X 


X 


-4.4* 


X 


X 


-13.1* 


-5.4* 


-3.3 


X 


X 


-14: 


X 


X 


2.3 


-5.7 


5.6* 


X 


X 


Jan 


AS 205 


-12.4 


X 


X 


-28.2 


X 


X 


-16.3* 


-7.5* 


-6.1 


-9.7 


375.1* 


-15.4 


-7.8 


X 






5.9* 


X 




AS 209 


11.4 


X 


X 


11.0 


X 


X 


-15.0* 


-5.2* 


-10.0 


4.6 


X 


-20.0 


-6.4 


X 


X 


25.7 


4.8* 


X 


X 


oo ■ 


AS 353 


-197.1* 


-28.3 


X 


-194.4* 


16.3 


X 


0.7 


-2.9* 


X 


0.3 


X 


X 


-10.6 


X 


X 


X 


X 


-21.4 


-14.4 


T— 1 1 


BP Tau 


62.9* 


X 


X 


53.7* 


X 


X 


-10.4* 


-1.6* 


-1.8 


0.4 


37.9* 


-2: 


25: 


X 


-7.7 


-9.0 


7.5* 


X 


X 


1 ■ 


CI Tau 


X 


X 


X 


X 


X 


X 


-16.0* 


-6.9* 


-3.6 


-1.0 


X 


-16.4 


-4.5 


X 






2.5* 


-27.5 




22v 


CO Ori 


-44.8* 


-13.4 * 


X 


-38.5* 


-5.2 * 


X 


-9.3* 


3.2* 


X 


X 


X 


X 


-66.0 


3.0* 


X 


X 


11:* 


X 


X 


CW Tau 


X 


X 


X 


X 


X 


X 


-5.7* 


1.6* 


X 


4.6 


X 


3.4 


13.0 


X 


-11.1* 


-3.7 


X 


-0.4 


X 


m ■ 
^ ■ 

o : 


DE Tau 
DF Tau 


X 

-65.2* 


X 

-13.1 


X 

141. f 


X 

) -62.2* 


X 

-23.2 


X 

136.4 


X 
X 


-2.8* 
4.1* 


6.9 
-1.4 


6.8 
-11.7 


X 
X 


-10.8 
-0.7 


1.5 
-7.2 


X 
X 


X 

-5.7 


3.4 
-1.6 


X 
X 


X 

-7.5 


2.4 
-5.5: 


ph/00 


DC Tau 


-95.4* 


-6.8 


X 


-91.9* 


6.9 


X 


-2.5 


-0.2* 


X 


4.9 


X 


X 


8.2 


X 


X 


10.2 


7.2* 


2.9 


5.9 


DK Tau 


-16.3* 


X 


X 


-11.9* 


X 


X 


-11.3* 


-2.2* 


-1.4 


93.7 


X 


-5.3 


91.2 


X 


200:* 


32.6 


6.7* 


X 


X 


DL Tau 


-122.7* 


-18.1 


X 


-128.1* 


13.2 


X 


-20.1* 


-6.2* 


-3.5 


-1.6 


X 


-6.0 


2.9 


X 


X 


15.1 


X 


-1.8 


X 


o ■ 


DN Tau 


-8.9 


X 


X 


1.9 


X 


X 


-10.9* 


-0.8* 


-2.7 


X 


X 


-9.1 


X 


X 


-5.5 


X 


8.0* 


X 


X 


astr 


DO Tau 
DQ Tau 


-84.6* 
-24.4 


-4.2 
100:* 


X 
X 


-80.0* 
-23.6 


1.0 

101:* 


X 
X 


X 

9.1* 


-0.2* 
4.0* 


-3.3 
0.3 


-6.4 
-45.6 


X 
X 


-5.1 
-1.4 


-3.6 
-54.9 


X 
X 


X 
X 


-3.1 

X 


X 
X 


-20.3 
-0.2 


X 
X 


v: 


DR Tau 


-9.74 


162.9* 


X 


-4.14 


108.7* 


X 


X 


-1.3* 


X 


-2.5 


350.7* 


X 


-4.4 


X 


X 


-3.1 


X 


-1.7 


-2: 




DS Tau 


-4.1* 


X 


X 


-4.0* 


X 


X 


-11.6* 


-7.9* 


-9.7 


X 


X 


-11.9 


X 


X 


-10.4 


-19.3 


7.2* 


X 


X 


:ar] 


GG Tau 
GK Tau 


X 
X 


X 
X 


X 
X 


X 
X 


X 
X 


X 
X 


-7.0* 

X 


1.0* 
-1.7* 


1.0 
-4.8 


50.7 
-30: 


X 
X 


-10: 
-5: 


4: 

X 


X 
X 


X 

-12.0 


X 

-42: 


8.2* 
3.5* 


X 
X 


X 
X 




CM Aur 


5.4* 


X 


X 


0.8* 


X 


X 


-11.0* 


-4.3* 


-3.0 


X 


X 


-6.7 


X 


X 


-4.1 


-15.4 


7.1* 


X 


X 




GW Ori 


-15.0 


X 


X 


4.9 


X 


X 


-2.8* 


3.4* 


X 


-1.6 


X 


X 


-1.6 


X 


X 


X 


12.0* 


X 


X 




HL Tau 


2.3 


X 


X 


6.3 


X 


X 


-8.6* 


-0.6* 


X 


-10.1 


X 


X 


-8.5 


X 






X 


-11.9 






RW Aur 


-93.0* 


12.1* 


X 


X 


-85.0* 


10.7* 


-13: 


0.2* 


X 


X 


X 


X 


-11.0 285.8* 


-52.0* 


37.3 


X 


X 


X 




RY Tau 


27.2* 


-98.0 


X 


40.4* 


-84.2 


X 


9.1* 


17.8* 


X 


-2.3 


X 


X 


5.3 


X 


X 


5.0 


27.1* 


X 


X 




SU Aur 


-4.3* 


X 


X 


0.2* 


X 


X 


-11.8* 


-3.4* 


X 


X 


11.5* 


X 


X 


10.6 


X 


X 


2.2* 


X 


X 




T Tau 


-161.6* 


-7.76 


X 


-161.9* 


-9.73 


X 


-14.6* 


-5.4* 


X 


-4.46 


X 


-13: 


-12: 


X 


X 


25.2 


5.7* 


-16.4 


X 




UX Tau 


-9.8* 


X 


X 


-5.9* 


X 


X 


-12.8* 


-3.2* 


-9.0 


X 


X 


-13.1 


X 


X 






6.3* 


X 






UY Aur 


X 


X 


X 


X 


X 


X 


X 


2.6* 


18.9 


-35.6 


X 


0.9 


-23.1 


X 


-10.6 


10.4 


X 


X 


X 




V410 Tau 


5.2* 


X 


X 


2.7* 


X 


X 


-1.3* 


7.8 


49.2 


X 


X 


46.5 


X 


X 


5.2 


X 


10.5 


X 


22.8 
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Table 1. Reflection velocity (km s ^) to match the red and blue wings. 





Star 


Ha 


H/3 


H7 


m 


Ca II 


Ca II 


o 
o 












A8662 


A8498 


o 
















an 


AA Tau 


-30 


ASb 


ASb: 


ASb 


Qa 




1—5 

OO 


AS 209 


2 


2 


10 





5 


-7 


T— 1 


AS 205 


0^ 








-10 


-10: 


> 


AS 353 





20: 


0: 


0: 





-10 


(N 
(N 
m 


BP Tau 


ASb^ 


ASb" 


ASb" 


ASb 








'0001 


CI Tau 


-10 











-5 


CO Ori 


23 


X 


X 


X 


X 


X 




CW Tau 





10 


X 


X 


5 


ASr'' : 


astro- 


DE Tau 
DF Tau 






12 
-18 


ASb" 
-27 


ASb: 


7: 
ASb: 




ASb: 


iv: 


DC Tau 


-33 


-33 


20 


20 


5 


ASr"^ 


1 


DK Tau 


ASr^ 


ASb: 


ASb 


ASb 


ASr: 


ASr: 




DL Tau 





8 


10 


20 





ASr"^ 




DN Tau 


AS* 


ASb^ 








Qa 


-8'^ 




DO Tau 





-3 


X 


X 


-6 


-4 




DQ Tau 


ASr'^ 


ASr"^ 


-6 


-10 


ASb 


ASb 




DR Tau 








-20 





ASr^ : 


ASr^ : 




DS Tau 


-35 


-50 


-30 


-35 




-12" 




GG Tau 


ASb^: 





5 


40 


ASr 
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Table 1 — Continued 



Star Ha Rp H7 US Call Call 

A8662 A8498 







00 






"4: . 




GM Aur 


ASb"^ 


-50: 


-30 


ASb 


0" 


-5" 


GW Ori 





-22 


X 


X 


ASb'': 





HL Tau 











ASb" 


-8 


RW Aur 


ASr^ 


10 


X 


X 


ASb" 


ASb 


RY Tau 


ASr'^' 


X 


X 


X 


ASr 


5 


SU Aur 


-13: 


X 


X 


X 


X 


X 


T Tau 





-5 


15: 





-5 


-10 


UXTau 


-42 








-8" 


-15" 


UY Aur 


ASb" 


-15 


-20 




ASb: 


ASb: 


V410 Tau 


ASb 


X 


X 


X 


X 


X 



Note. — Lines identified as "■ • ■" were not available 
and as "x" were not reliably measured or observed in 
absorption. ASb means that the line is blue ward asym- 
metric and ASr redward asymmetric. "=only the nar- 
row component is present, ''=only the low velocity part 
is symmetric, the wings at high velocities are asymmet- 
ric, "=only the high velocity is symmetric, the rest of 
the line is asymmetric, '^=the high and low velocities are 
asymmetric, the line part between those is symmetric. 



Profiles of Strong Permitted Lines in Classical T Tauri Stars 



Silvia H. P. Alencar^'^ and Gibor Basri^ 



ABSTRACT 

We present a spectral analysis of 30 T Tauri stars observed with the Hamilton echelle spectro- 
graph over more than a decade. One goal is to test magnetosphcric accretion model predictions. 
Observational evidence previously published supporting the model, such as emission line asym- 
metry and a high frequency of redshifted absorption components, are considered. We also discuss 
the relation between different line forming regions and search for good accretion rate indicators. 

In this work we confirm several important points of the models, such as the correlation 
between accretion and outflow, broad emission components that are mostly central or slightly 
blueshifted and only the occasional presence of redshifted absorption. We also show, however, 
that the broad emission components supposedly formed in the magnetospheric accretion flow 
only partially support the models. Unlike the predictions, they are sometimes redshifted, and 
are mostly found to be symmetric. The published theoretical profiles do not have a strong 
resemblance to our observed ones. We emphasize the need for accretion models to include a 
strong turbulent component before their profiles will match the observations. The effects of 
rotation, and the outfiow components, will also be needed to complete the picture. 

Subject headings: line: profiles — stars; pre-main-sequence — stars: formation 



1. Introduction 

Classical T Tauri stars (CTTSs) are young, al- 
most solar-mass stars that exhibit a wide range of 
permitted and sometimes also forbidden emission 
lines, together with an excess continuum emis- 
sion that goes from the infrared to the ultravio- 
let. Their spectral energy distribution is consis- 
tent with the presence of a circumstellar disk that 
appears to play a major role in the regulation of 
both the infall and outfiow of material from the 
star-disk system. A strong wind component is also 
thought to be present and cause the blueshifted 
absorption features commonly seen in the Balmer 
lines. 

In the past years many attempts have been 
made to explain those characteristics. The per- 
mitted line emissions were discussed in terms of 
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outfiowing winds (Hartmann et al. 1982, 1990; 
Natta & Giovanardi 1990), turbulence in a bound- 
ary layer between the disk and the stellar surface 
(Bertout et al. 1988; Basri & Bertout 1989) 
and chromospheric activity (Calvet et al. 1984, 
1985). More recently, magnetospheric accretion 
models were proposed (Hartmann et al. 1994; 
Shu et al. 1994; Hartmann 1998). In these mod- 
els the circumstellar disk is truncated by stellar 
magnetic field lines and the strong broad emis- 
sions arise from the accelerated infalling material 
channeled through the lines that connect the disk 
to the star. At the base of the magnetospheric ac- 
cretion column, where the accreted material hits 
the star, hot spots and the hot continuum emission 
(veiling) are produced. 

Edwards et al. (1994) showed observational ev- 
idence in a sample of 15 CTTSs that they claimed 
confirm some of the magnetospheric predictions, 
like blueward asymmetric emission lines (due to 
disk occultation of the receding fiow) and red- 
shifted absorption components at typical free-fall 
velocities (inverse P Cygni profiles), which are nat- 
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urally explained by the magnetospheric infall. The 
line profiles of several stars were also well repro- 
duced by magnetospheric calculations (Muzerolle 
et al. 1998a; Najita et al. 1996), but often the 
theoretical profiles are more asymmetric and al- 
ways less broadened than the observed ones. How- 
ever, wc must keep in mind that the models of 
Hartmann et al. (1994) do not include rotation 
and winds that would certainly influence both the 
asymmetry and the broadening of the theoretical 
profiles. 

The CTTSs were initially found to be slow ro- 
tators due to observational surveys that reported 
a bimodal distribution of stellar rotation periods 
among TTS (Choi & Herbst 1996), where the 
slow rotators (P > 4 days) seemed to exhibit disk 
signatures, as near infrared excess emission, that 
the fast rotators did not (Edwards et al. 1993). 
It has always been a challenge to imdcrstand the 
spin-down of those stars but the magnetospheric 
model predicts that the magnetic interaction be- 
tween the star and the disk could regulate the stel- 
lar rotation. This could explain why weak T Tauri 
stars (WTTS), that do not have disks present a 
wider range of rotational periods. Recently Stas- 
sun et al. (1999), analyzing the rotation period 
distribution of young stars in the Orion nebula re- 
gion, showed that this bimodal distribution does 
not seem to exist and that the earlier works con- 
clusion was also statistically consistent with a uni- 
form distribution. Herbst et al. (1999) suggest 
this recent result, however, may be confined to 
low mass TTSs. 

In view of the many points still under discus- 
sion, we present in this paper the spectral analy- 
sis of an observed sample of 30 TTSs with a large 
spectral coverage that are used to analyse the mag- 
netospheric model's predictions and test the pre- 
vious observational results. The observations and 
reduction procedures are described in Section 2, 
the equivalent width and veiling measurements are 
detailed in Section 3, we analyse the many features 
in our spectra in Section 4, discuss the results in 
Section 5 and draw our final conclusions in Section 
6. 

2. Observations 

We present the spectral analysis of a sample of 
30 TTSs listed in Table 1. The observations, which 



span over more than a decade, were carried out at 
Lick Observatory with the 3m Shane telescope and 
the Hamilton Echelle Spectrograph (Vogt 1987) 
coupled either to a TI 800x800 CCD or a FORD 
2048x2048 CCD. With the smaller detector two 
settings were established: a red setting covering 52 
orders from 4900A up to 8900A and a blue setting 
covering 38 orders from 3900 A up to 5200A. The 
spectral coverage, however, is not complete, with a 
gap of a few dozens of angstroms between orders. 
Whenever possible, blue and red observations were 
obtained in the same night or within 1 night of 
difference, but for some stars wc only have the red 
setting. The bigger CCD installed in 1992 permits 
a full spectral coverage of ~ 92 orders ranging from 
3900 A to 8900 A. The mean resolution of the 
spectrograph is A/AA ~ 48, 000 and the exposure 
times varied from 15min to IhlSmin, depending 
on the target and on the CCD used. 

The reduction was performed in a standard way 
described by Valenti (1994) which includes flat- 
fielding with an incandescent lamp exposure, back- 
ground subtraction and removal of cosmic rays. 
Wavelength calibration is made by observing a 
Thorium- Argon comparison lamp and performing 
a 2-D solution to the Thorium lines. Radial and 
barycentric velocity corrections are applied and all 
the data shown are in the stellar rest frame. The 
spectra are not flux calibrated, so each spectrum 
has been continuum normalized. Due to differ- 
ences in weather conditions, exposure times and 
efficiency between different chips there is a wide 
range of signal to noise in the data. However, not 
only strong emission line profiles were reliably ex- 
tracted but also many absorption line profiles as 
shown in Figure 1. After the reduction procedure 
the spectra have been binned and smoothed with 
a window width of 3. A weak-line T Tauri star 
was included in the sample, V410 Tau (Herbig & 
Bell 1988), as a reference of stars without disk ac- 
cretion. DQ Tau is a known binary system (Basri 
et al. 1997), the chosen spectrum corresponds to 
an outburst event when the lines are stronger and 
to a phase (</> = 0.02) where there is no separation 
between the lines of the two components. 
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3. Measurements 

3.1. Equivalent widths and profile decom- 
position 

Wc first measure the equivalent widtli of all the 
lines in Figure 1 that could be reliably identified. 
The results are presented in Table 1 and will be 
used in the following sections to study the rela- 
tions between "fluxes" (veiling corrected equiv- 
alent widths) from different lines and between 
fluxes and parameters such as the mass accretion 
rate. 

The line profiles were also decomposed with a 

multiple gaussian procedure based on the Mar- 
quardt method, so that we could also study cor- 
relations between line components. Most of the 
lines were very well fitted by 1 or 2 gaussians 
and only sometimes required a third component. 
Some exceptions are the flat topped profiles of T 
Tau, SU Aur and DE Tau that are clearly non- 
gaussians. The fitting procedure was not always 
straightforward, since some lines are contaminated 
by non-stellar emission, some are blended by other 
lines and some can be fitted with more than one 
gaussian combination (like 2 emissions instead of 
1 emission and 1 absorption, for example). We 
explain below how we dealt with those problems. 

The Ha line often exhibits a double peaked pro- 
file which was, in most cases, better fitted by a 
broad emission and a blue absorption than 2 emis- 
sions. The Ha wings are not steep and seem to 
belong to a very broad emission line rather than 
to 2 narrow ones. It is often possible to confirm 
that a true absorption is present by looking at the 
upper Balmcr lines where the feature is enhanced 
relatively to the emission and dips below the con- 
tinuum. 

The Na D lines are strongly contaminated by 
the vapor lamps of the city of San Jose but 
our decomposition procedure was able to reliably 
subtract that contamination with a simultaneous 
gaussian fit while decomposing the line profiles. 

The Ca II and He I lines exhibit profiles typi- 
cally with narrow and broad emission components; 
some lines show only one of the components. In 
some Ca II lines the narrow component is present 
inside a photospheric absorption core and, in that 
case, the emission-line equivalent width was mea- 
sured from the base of the absorption core in- 



stead of the continuum. The Ca II infrared triplet 
(IRT) lines are also blended with a Paschen emis- 
sion line that has usually strengths less than 20% 
of the IRT, as also noted by Hamann & Persson 
(1992). Whenever the Ca II lines are weak or nar- 
row, the Paschen contamination is significant and 
has been extracted from the equivalent width mea- 
surements. 

Finally, the Fe II equivalent widths were only 
measured when the line was in strong emission 
and could surely be distinguished from the lines 
nearby. 

3.2. Veiling 

The equivalent widths presented in Table 1 are 
contaminated by the variable continuum against 
which they were measured. They need to be cor- 
rected for the veiling in order to represent a true 
emission strength that can be regarded as a "flux" 
measurement. Some of our spectra have measured 
veilings in previous papers (Basri & Batalha 1990; 
Johns-Krull & Basri 1997) but in order to be 
consistent we decidcid to measure; all the; veilings 
again. The new values are presented in the last 
columns of Table 1 and the common measurements 
agree within the errors with the previously pub- 
lished ones. 

In order to calculate the veilings we chose spec- 
tral orders with many photospheric absorption 
lines both in the blue and red settings. The red or- 
ders correspond to the wavelength ranges 5557A< 
A < 5642A and 5668A< A < 5755A and the blue 
orders to 457lA< A < 464lA and 4763A< A < 
4836A. The idea is to compare the spectrum of 
a CTTS with that of a standard star with the 
same spectral type that has been broadened to the 
CTTS rotational velocity and veiled. The choice 
of suitable standard stars is important and we fol- 
lowed, whenever possible, the suggestion of Basri 
& Batalha (1990) to use the Hyades dwarfs as re- 
liable standards. We also used the standard and 
CTTS spectral types determined by them. 

We first rotationally broadened the continuum 
normalized spectrum of the standard star and then 
applied the veiling to the standard's spectrum, 
comparing it with the observed CTTS spectrum 
until a good match was found. The best adjust- 
ment was determined by eye, but a unique veil- 
ing value did not perfectly match all the photo- 
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spheric linos and the results presented in Table 1 
correspond to a mean value obtained with all the 
lines in the selected orders. The errors in the veil- 
ing determination arc bigger for the higher veiling 
values, as the photospheric lines of highly veiled 
stars are extremely shallow and hard to match. 
Typically wc have errors of 0.1 for veilings smaller 
than 1.0 and up to 0.5 for the higher ones. Most 
of the rotational velocities were taken from Basri 
& Batalha (1990), although we verified them by 
checking the match between the FWHM of the 
cross-correlation normalized profiles of the origi- 
nal standard with the CTTS spectra and of the 
original standard with the broadened and veiled 
standard. 

The veiling values obtained are used to search 
for correlations between veiling corrected line 
equivalent widths, defined as in Johns-KruU & 
Basri (1997): 

= W,^{V + 1) (1) 

where Weq is the measured equivalent width and 
V is the veiling. 

4. Analysis 

4.1. General Characteristics 

The magnetospheric accretion model is the cur- 
rent consensus model to describe the accretion 
processes in CTTSs. These models predict gen- 
eral trends like central or blueshifted asymmetric 
broad emission lines, sometimes with redshifted 
absorption components (Hartmann et al. 1994; 
MuzeroUe et al. 1998b), that we should be able 
to verify in our large sample of stars and lines. 
The magnetospheric accretion scenario may also 
include winds that would cause further emission 
and blueshifted absorption (Shu et al. 1994). 
Some theoretical predictions have been confirmed 
by observational results (Edwards et al. 1994; 
Muzerolle et al. 1998a) but there is still a long 
way to go in order to fully understand the nature 
of the processes taking place in CTTSs. Below are 
some general profile characteristics from our sam- 
ple of stars that we will analyse in detail in the 
next sections. 

Many of our emission lines show broad com- 
ponents (BC) that tend to have blueshifted cen- 
troids (Ha, H/3, the He I BC, the IRT BC, the 



Na D and the Fc II A4923) following the magne- 
tospheric predictions, but the centroids are also 
redshifted in some stars and exhibit a wide range 
of shift values. The broad emission components 
are also generally well fitted by gaussians which 
may have a turbulent origin (Basri 1990) rather 
than the infall scenario. Redshifted absorption is 
not commonly present. The clearest example is 
BP Tau, where it appears in the Balmer lines, the 
Na D and the Ca II (A8662). Blueshifted absorp- 
tion is a much more common feature that can be 
seen in many Balmer and Na D lines and some- 
times in lines like Ca K in RW Aur. Many lines 
show a multiple-component profile that can of- 
ten be decomposed either as narrow and broad 
components, e.g. He I, Ca and Fe lines (Basri & 
Batalha 1990; Beristain et al. 1998), or as emis- 
sion and absorption components, like the Balmer 
lines. The different components often show dis- 
tinct characteristics (like broadening and shifts) 
suggesting that the lines may originate in distinct 
regions. The Fe I lines are mostly in absorption 
but can be found in strong emission in very active 
stars (as in AS 353, DG Tau and RW Aur) always 
with a double-peaked profile that could bo due to 
NLTE chromospheric core inversion (inversion of 
the source function) or disk rotation. 

4.2. Line Strengths 

The study of the correlations between the fluxes 
of different lines is an important tool that gives 
insight on the relations between different line for- 
mation regions. Correlations can indicate common 
forming regions but can also mean that a common 
process is affecting both regions at the same time. 
In this section we discuss the correlations that ap- 
pear and also those that do not in our data. 

The Balmer lines emission components arc well 
correlated with each other, as are the IRT lines 
with each other, which is expected as they should 
be formed in regions with the same physical char- 
acteristics. The difference between them should 
arise mainly due to differences in line optical 
depths. In the Balmer lines case, the correlation 
seems to decrease as we compare lines that are fur- 
ther apart in the series, e.g. Ha and H/3 show a 
strong correlation with each other (Figure 2, top) 
with a false alarm probability (FAP) of 3.24 x 10"^, 
but Ha and Rd do not (FAP=0.168). This behav- 
ior is biased by stars with strong wind contribu- 
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tion, where the emission component is sometimes 
almost absent due to large absorption components 
and/or a hot underlying photosphere. Eliminat- 
ing those stars the good correlations are partially 
recovered (Figure 2, bottom), the FAP of the 
Ha and RS emission components decreasing to 
3.55 X 10"'' and presenting a linear correlation co- 
efficient r=0.880. 
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Fig. 2. — Veiling corrected equivalent widths of 
the emission components: vs. W^g and vs. 



The entire Ha and He I lines are also well corre- 
lated (Figure 3 ; FAP=1.26 x IQ-^), except for DQ 
Tau which has a very high Ha equivalent width in 
the outburst spectra we chose and CW Tau which 
has a high uncertainty in the veiling measurement 
(we get FAP=4.65 x 10-i° without DQ Tau and 
CW Tau) . The same kind of correlation is verified 
with H/3, but H7 and US do not show a well defined 
correlation with He I. The apparent lack of correla- 
tion with H7 and H^ may be due to the small num- 
ber of spectra where those lines could be reliably 
measured and the lower S/N present in our spec- 
tra in that region. Apart from the correlation with 
the entire lines, we also verified some good correla- 
tion between the Ha emission component and the 
He I BC (FAP=4.68 x IQ-^ with DQ Tau and CW 
Tau and FAP=8.32 x lO""' without them). The 



He I line is believed to form in high temperature 
regions (20,000K to 50,000K) and together with 
the large width of its BC it suggests that the BC 
comes either from the highly turbulent shock re- 
gion or the magnetospheric flow, where the broad 
emission component of the Balmer lines is also 
supposed to be formed. Due to their larger op- 
tical depths, the Balmer lines probe a much larger 
volume of gas than the He I lines, but the correla- 
tions between the entire lines and the components 
show that their formation regions are strongly re- 
lated to each other. 

Also thought to be formed in the accretion flow 
because of its large widths and blueshifted cen- 
ter, the broad IRT component does not strongly 
correlate with the He I BC (FAP=1.97 x 10"^ 
and r=0.678 for A8498) or with the Ha emission 
(FAP=6.12x 10-3 and r=0.590 for A8498). For the 
He I BC, no correlation is found for small equiva- 
lent width values and a tendency for the equivalent 
widths to increase together exists, although with 
much scatter. For the Ha, however, a correlation 
seems to exist except for the higher Ha flux values 
(where the scatter is really large). Although the 
Ca II line is formed at much lower temperatures, 
it is expected to strongly correlate with Ha and 
He I if the accretion flow is sufficiently organized 
that a global change in accretion rate is felt by all 
diagnostics. 

Batalha et al. (1996) showed that the nar- 
row components (NC) of CTTSs have a flux ex- 
cess when compared to the WTTSs and this was 
attributed to the reprocessing of radiation pro- 
duced in accretion shocks as the accreting ma- 
terial hits the stellar atmosphere. However, nei- 
ther the He I NC nor the Ca II NC correlate with 
the Balmer emission, which implies that although 
both the NCs and the Balmer lines are related to 
accretion, they are probably probing regions that 
are far away from each other and consequently not 
equally influenced by the same physical processes, 
or the time delay for changes that affect one region 
to reach the other is too large for the changes to 
show up in simultaneous spectra of both compo- 
nents. 

We also confirmed the trend noticed by Batalha 
et al. (1996) that the IRT and the He I NCs are 
positively correlated (Figure 4 and FAP=1.41 x 
10"^), (unsurprisingly, since we have a lot of data 
in common), while MuzeroUe et al. (1998a) found 
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that they were uncorrelated. Those hnes are gen- 
erally thought to sample quite different temper- 
ature regions - the He I is formed in high tem- 
peratures (> 20, OOOK) while the IRT is a chro- 
mospheric line formed at 4, OOOK < T < 7, OOOK. 
The correlation foimd is probably related to the 
excess emission due to the influence of the accre- 
tion shock on the stellar atmosphere, a common 
feature in both of them. 

Iron lines are quite common in CTTSs, so we 
decided to analyse the behavior of four of them (2 
Fe I and 2 Fe II) in our sample. The Fc II (A4923) 
line was only fitted when in emission. Only 40% 
of the sample presented this line in emission and 
we were able to measure the veiling for just some 
of them, but before applying the veiling correc- 
tions it shows a good correlation with the He I BC 
and with the IRT BC. These correlations may in- 
dicate that the Fe II lines, at least when in strong 
emission, share a common emitting region with 
the BCs of He I and the IRT. 

The Fe I absorption lines (A6192 and A5497), 
with excitation energies of 4.43eV and 3.27eV re- 
spectively, are not correlated with each other. Al- 
though visually they both have some resemblance 
to the Li I line, nothing was found between the 



Fe I (A6192) and Li I, while a strong correlation 
exists between the Fe I (A5497) and Li I lines. 
This is due to a temperature effect: the cooler 

stars show the deepest lines. 

Another important aspect of the correlation 
measurements is the quest for good accretion rate 
indicators that would later provide an easy way 
to obtain these rates for other stars. Using the 
mass accretion rates in Hartmann et al. (1998), 
there is a very good correlation between the IRT 
veiling corrected equivalent widths (Ca II A8498 
and A8662) and mass accretion rates (Figure 5 
(top) and FAP=1.80 x 10"''), as in MuzeroUe et al. 
(1998a) for Ca II (A8542). The only exception is 
DQ Tau, that has a log M = -9.4 MQyr"^ accord- 
ing to Hartmann et al. (1998) but —7.3 M0yr~^ 
according to Hartigan et al. (1995) (the largest 
discrepancies between their values). DQ Tau is 
known to have a variable accretion rate that de- 
pends on the orbital phase of the binary system 
and on events such as outbursts that may occur 
when the stars approach each other (Basri et al. 
1997). We decided not to use this star. 

The Balmer lines and their components are also 
correlated to the mass accretion rates (Figure 5 
(middle and bottom) with FAP=3.40 x lO"'' and 
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7.09 X 10~^ respectively). The results for the ab- 
sorption component are not as reliable as the oth- 
ers, mainly due to the small number of stars that 
exhibit a clear absorption and for which there are 
published mass accretion rates. The flux of the 
emission and absorption components increases as 
the accretion rate increases, showing that the ac- 
cretion may be powering both the emission and the 
outflow. This would be consistent with theoreti- 
cal results (Shu et al. 1994) that suggest that the 
wind mass loss rate and the disk accretion rate are 
directly proportional (but see below for caveats). 

4.3. Evidence for Outflows 

Winds are expected to be present in CTTS and 
help carry away angular momentum. If cool, they 
should appear mainly as absorption components. 
Due to the presence of the disk the absorption is 
expected to be blueshifted, the receding part of the 
wind is blocked by the disk from the observer's line 
of sight. The emission components of CTTS were 



thought for some time to be produced in optically 
thick winds, but these models also generated up- 
per Balmer line profiles with very deep absorptions 
that were not consistent with observations. 

The Balmer and Na D blueshifted absorption 
components are the most clear evidence in Figure 
1 that a strong wind is present in our sample of 
CTTSs. Almost 80% of our stars show blueshifted 
absorption components in at least one line, the 
most common being Ho; (see Table 2). The num- 
ber of blueshifted components decrease as we look 
at the upper Balmer lines, appearing only in ~ 
35% of the profiles (where they tend to be 
weak). This suggests that the blueshifted absorp- 
tions arise from a region that is optically thin at 
least in the upper Balmer lines, and distinct from 
the emission component that is believed to arise 
mainly in the optically thick accretion flow. Some 
stars, like DR Tau and AS 353, present P Cygni 
profiles in all the Balmer lines. These are normally 
associated with strong winds, and other stars, like 
CO Ori, RW Aur, SU Aur and RY Tau, have such 
a strong absorption component that it sometimes 
suppresses the emission in the upper Balmer lines 
of the stars with the hottest underlying photo- 
spheres. A number of these stars are known to 
be jet sources. 

The 8 stars in our sample that exhibit blueshifted 
absorption components in the Na D lines (presum- 
ably indicating high mass loss rates; Hartmann 
1998), also do in Ha. The lack of P Cygni pro- 
file in the Na D lines does not necessarily mean 
low mass loss rates, as the star may present a 
deep photospheric absorption that prevents the 
detection of the blue absorption due to the wind 
(Natta & Giovanardi 1990). Due to the difference 
in their optical depths, the Balmer and Na lines 
will sample different wind volumes. As resonance 
lines, the sodium probes the coolest part of the 
wind. The blueshifted components of both lines 
tend to be centered at approximately the same 
velocities (Table 2), likely meaning that different 
wind regions share some common kinematic com- 
ponent. The Na D blueshifted absorptions bear 
no other resemblance to the Ha ones, being much 
weaker in intensity and width. The position of 
blueshifted absorption centers in these lines vary 
from star to star suggesting that the wind veloc- 
ities can occur with a wide range of values (from 
to -270 km s~^ in our sample, see Figure 6). 
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Fig. 6. — Blueshifted absorption component's 
peak position. 



There may well be projection effects that must be 
taken into account; these would imply the outflow 
is not spherical. 

After decomposing the profiles we can also in- 
vestigate the correlations that appear between dif- 
ferent wind components and between wind com- 
ponents and other features in order to better un- 
derstand the wind forming region. The Ha and 
H/3 blueshifted absorptions are very well corre- 
lated (FAP=2.39 X 10"^), showing that they both 
probe the same wind region, and although there 
is apparently no correlation with the blueshifted 
absorptions of H7 and H^, this may be due to the 
small number of cases. Where the absorption does 
appear in the upper lines the wind component is 
often very weak and consequently hard to fit. 

The Ha and H/3 blueshifted absorptions also 
show some correlation with their line's emis- 
sion components: the stronger the emission, the 
stronger the absorption. Unfortunately, that does 
not have a unique interpretation. One possibility 
is that the same mechanism that is powering the 
emission is also powering the wind in the sense 
that when more material is accreted there is more 
to outflow. This is predicted by the accretion 



model and also means that the absorption region 
is outside the emission region. A second possibil- 
ity is that if both emission and absorption were 
formed in winds, the absorption component would 
be due only to the decay of the source function 
in the outer regions of the atmosphere and not to 
the accretion itself. If the first possibility is cor- 
rect, this relation could also explain the trend seen 
between the Ha emission and absorption compo- 
nents and the mass accretion rates. 

Finally, the present result may also be an effect 
of line optical depth; if the wind is optically thick 
the central depth of the blueshifted absorption 
is fixed and the changes that appear are in fact 
only due to the emission component variations as 
shown by Johns & Basri (1995b) in the case of SU 
Aur. Johns-Krull & Basri (1997) also found that 
the Ha absorption and emission strengths of many 
DF Tau spectra were correlated and that the ab- 
sorption component strength correlated with the 
veiling. 

A curious fact about the absorption compo- 
nents can be noticed in the stars that show an 
absorption component in the 4 analysed Balmer 
lines. This component tends to be shifted towards 
the red as we go from Ha to US and in Figure 
7 we plot the absorption centers vs. the Balmer 
lines oscillator strengths for those stars. It seems 
that we may be seeing the acceleration of the flow. 
This shift of the absorption component can also 
be noticed in Edwards et al. (1994) in the resid- 
mil profiles of DK Tau and in Appenzeller et al. 
(1988) in the DR Tau spectra, showing that this 
behavior may be quite common among CTTSs. 

Another feature that is associated with strong 
activity, presumably both inflow and outflow, are 
the iron lines in strong emission. These lines are 
often in absorption in our sample but Fe II (A4923) 
is also found in emission in ~ 40% of our stars. 
When this happens, a wind component is present 
at least in Ha, sometimes also in the other Balmer 
lines and in the Na D lines, and the IRT lines al- 
ways show a BC. The Fe lines in emission seem 
to be related to both the presence of outflow and 
high veiling, as some stars that show blueshifted 
absorption in Ha but no BC in the IRT do not 
present Fe II (A4923) in emission. DG Tau, AS 353 
and RW Aur, that also exhibit Fe I (A6192) and 
Fe II (A4352) in strong emission, also have associ- 
ated jets (Edwards et al. 1994; Mundt & Eisloffel 
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Fig. 7. — Balmer line absorption center vs. os- 
cillator strength. RW Aur, DR Tau and DF Tau 
are plotted twice as 2 wind components were mea- 
sured. 



1998). But unlike the forbidden lines that tend to 
be blueshiftcd, the iron lines when in emission are 
centered at the stellar rest frame. They are prob- 
ably not produced in the outflow (wind or jet), 
since the occultation of the receding part of the 
flow by the disk would yield a blueshifted line. 

4.4. Evidence for Infall 

Redshifted absorption components in the line 
profiles (inverse P Cygni profiles), can be the clear- 
est spectral evidence that infall of material is oc- 
curring, and are naturally predicted by the mag- 
netospheric accretion models. These models, how- 
ever, require special conditions for the redshifted 
absorption component to bo visible, such as low 
line thermalisation, low emission damping wings, 
and proper inclination values. The inclination de- 
pendence arises from the contrast of the line source 
function of the infalling gas and the continuum 
source function where the gas stream is projected 
(Hartmann et al. 1994). If the projection is 
against the cool photosphere (low inclinations) no 
absorption is found, but if it is against the hot 
shock material, like a ring or spots, (higher incli- 



nations) the absorption can be produced. 

Edwards et al. (1994) presented an analysis 
of 15 CTTSs included in our sample, and showed 
that 87% of their stars had redshifted absorption 
components in at least one line at ~ 200 to 300 
km s~^, a range of velocities that is consistent with 
the ballistic infall predicted by the magnetospheric 
models for a typical CTTS. Only ~ 40% of our 
stars have redshifted absorption components in at 
least one line and generally not at ~ 200 to 300 
km s~^. BP Tau, for example, shows redshifted 
absorption in many lines at ~ 50 km s~^ and UX 
Tau has a strong redshifted absorption in Ha at 43 
km s~^. The difference in the absorption centroid 
position of the stars from the values predicted by 
the models could partially be due to inclination 
effects. The difference in the number of redshifted 
absorption occurrences found in the two different 
analysis could be due to the fact that Edwards et 
al. (1994) analysed residual profiles instead of 
the original ones. This gives, according to them, 
a high sensitivity in defining features such as the 
Balmer lines wings, where the redshifted absorp- 
tion lies. However, in the residual profiles showed 
by Muzerolle et al. (1998a) of 11 CTTSs, red- 
shifted absorption at typical free-fall velocities can 
be found in less than half of their sample, which 
is more in agreement with the predictions of the 
magnetospheric accretion models. 

In order to compare our results with the ones 
cited above we decided to generate residual pro- 
files from our sample, choosing 8 lines commonly 
found in emission: Ha, NaD, He I, Call (A8498 
and A8662), H/3, H7 and Fell (A4352). m was 
excluded as most of our spectra present low signal 
to noise in that region, making it hard to obtain 
reliable information about small absorption fea- 
tures. We generated the residual profiles by sub- 
tracting from each observation a broadened and 
veiled standard spectrum, using as standard stars 
the ones previously selected for the veiling calcu- 
lations. 

Although the residual method often enhances 
features like weak emission lines (sec Figure 8, 
top), care must be taken during the subtraction, 
as a small difference in continuum normalization 
between the two spectra can also introduce false 
features. The uncertainty in the veiling determi- 
nation and the diflaculty of choosing suitable stan- 
dard stars are other factors that can introduce er- 
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rors in the process. Finally, the legitimacy of the 
concept itself is dubious when the overlying emit- 
ting region is optically thick over most of the pro- 
file. 
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Fig. 8. — Residual profile determination, the thin 
solid line corresponds to the original spectrum, 
the dotted line to the standard spectrum and the 
thick solid line to the residual CTTS profile. Top: 
Fc II (A4352) is dramatically enhanced by the 
residual method. Bottom: the redshifted absorp- 
tion that appears in the unsubtracted spectrum 
in the Na D lines turns out to be mainly photo- 
spheric. 

The residual profile calculation suppresses, in 
principle, all the photospheric contribution that is 
present in the spectra, generating smoother pro- 
files, as the broad lines are no longer superimposed 
with narrow absorption ones. However, we no- 
ticed that it did not a priori enhance the redshifted 
absorptions present in our unsubtracted spectra; 
indeed some of the redshifted absorption compo- 
nents turned out to be mainly photospheric (for 
example, the Na D lines of BP Tau, DF Tau, DQ 
Tau and RY Tau, see Figure 8, bottom). In gen- 
eral, most of the redshifted absorptions became 
shallower, although less noisy, and we did not find 
many new ones. Our residual profiles show that 

50% of our stars present at least one line with 



a redshifted absorption component, a result in 
agreement with the original (unsubtracted) data, 
with the magnetospheric models, with the results 
by MuzeroUe et al. (1998a) but not with Ed- 
wards et al. (1994). Although very weak features 
such as the iron emission lines and the redshifted 
absorptions at free-fall velocities arc strongly af- 
fected by the underlying photospheric spectrum, 
one may notice that the general characteristics of 
the line profiles discussed in this work are weakly 
affected by the subtraction of the photospheric fea- 
tures and that our results are also valid for residual 
profiles. 

Finally, in agreement with both the model pre- 
dictions and with Edwards et al. (1994) we see 
less evidence of redshifted absorption in Ha at 
high velocities than in the upper Balmer lines, 
probably due to thermalisation effects and the 
Stark damping of the Ha wings. 

4.5. Profile Symmetry 

The magnetospheric accretion models predict 
that the EC emission is formed in the magneto- 
spheric infalling material. It should then exhibit a 
central or slightly blueshifted peak and blueward 
asymmetry due to oecultation by the disk and star 
of part of the flow moving slowly away from the 
observer. By blueward asymmetry we mean that 
there is an excess of blue emission compared to the 
red emission. Some theoretical profiles published 
by Hartmann et al. (1994) also show redshifted 
absorption at free-fall velocities, while others show 
extended red wings caused by the occultation of 
the high velocity blueshifted material by the star; 
both effects add up to enhance the line asymmetry. 

One way to search for asymmetries is to reflect 
one of the line wings about the broad emission cen- 
troid, as done by Johns-KruU & Basri (1997) for 
DF Tau. If the wings match each other, the line is 
considered symmetric, and this avoids a confusion 
between lines that are only shifted and lines that 
are really asymmetric. P Cygni profiles are more 
difficult to evaluate, since the blue emission wing 
is sometimes completely absent and it may be hard 
to determine the correct shift to be applied. 

In Table 3 we show the velocities that produced 
the best agreement between the opposite wings, 

and we also indicate which parts of the line are 
symmetric and which are not in the case of par- 
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tially symmetric profiles. The BCs that showed 
blueward asymmetry are labelled "ASb" in Table 
3, but they comprise only ~ 20% of the analysed 
lines. Most of the BCs are symmetric (well fitted 
by a Gaussian) and ~ 75% of them have a central 
or slightly blueshifted peak. In Figure 9 we show 
the velocity position of the broad emission peak for 
several lines and we notice that, although most of 
them are indeed blueshifted, a substantial num- 
ber present a rcdshifted peak. This is not seen in 
the published theoretical magnetosphere profiles. 
Lines with P Cygni profiles and those where the 
emission component was not reliably extracted are 
not included in this plot. 

Among the lines that presented asymmetries, 
the ones identified as "ASb*"' have symmetric low 
velocity material and a blueward excess of high ve- 
locity emission, and those identified as "ASr" show 
redward asymmetries. Neither type is predicted 
by the magnetospheric model. Some lines do 
present profiles with the characteristics described 
by the magnetospheric models, corresponding to 
"ASb" and "ASb'^" in Table 3. Examples of pro- 
files in agreement and opposite to the magneto- 
spheric symmetry predictions are shown in Figure 
10 with the reflected wings overplotted. 

Each Ca II IRT line is blended with a Paschen 
emission line and this could lead to a wrong asym- 
metry classification. However, from the thirteen 
stars that presented IRT asymmetries only four 
could be due to the Paschen contribution, which 
affects the high velocity wing redward of the line 
center. These stars are quoted as uncertain asym- 
metry measurements in Table 3. 

In general the symmetry predictions of the mag- 
netospheric model cannot be strongly confirmed 
with our sample. Contrary to them, most of the 
BCs are found to be symmetric; among the asym- 
metric ones many are not blueward asymmetric, 
and ~ 20% of the broad emission components 
present centroids redshifted by more than 5km 



5. Discussion 

Many aspects of the line formation of a CTTS 
can be inferred from the correlations found with 
the veiling corrected equivalent widths in the pre- 
vious sections. 

We showed that lines coming from a common 
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Fig. 9. — Broad emission component's peak posi- 
tion. Lines with P Cygni profiles and those not 
reliably measured are not included. 



lower level arc very well correlated, like II/3 with 
Ha and the IRT lines with each other. This is 
expected, as they need the same physical condi- 
tions to be produced and tend to be affected by 
the same physical processes. As previous stud- 
ies also showed, most lines have more than one 
component (BC and NC for the He I and the 
Ca II and emission and absorption components 
for the Balmer and Na D lines). These are likely 
formed in distinct regions, as the characteristics 
of those components are very different from each 
other throughout our sample and many correla- 
tions found were valid for one of the components 
but not for the other. According to our data the 
Balmer emission, the He I BC, the IRT BC and 
the iron lines (when in strong emission) are mainly 
formed in the magnetospheric flow. Alternatively, 
the He I BC may also come from the shock region, 
where the temperatiires are high and more appro- 
priate for the formation of that line instead of the 
magnetospheric flow itself. The NCs of the Ca and 
He I lines are thought to be formed at the stellar 
surface where it is perturbed by the infall of mate- 
rial and partly in the normal stellar chromosphere. 
The wind is the origin of the blueshifted absorp- 
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Fig. 10. — Emission lines with the red wings (dot- 
ted lines) reflected. The top lines arc symmetric 
and the bottom ones are asymmetric and in agree- 
ment with the magnetospheric model predictions. 



tion commonly seen in the Na D and Balmer lines. 

We also verified that the effect of accretion on 
the line components sometimes probably produces 
the correlations, as lines formed in regions with 
physically distinct properties (very different tem- 
peratures, for example) but both affected by ac- 
cretion can be found to correlate very well, like 
the He I and the IRT NCs. 

Somewhat puzzling is the lack of strong corre- 
lation between the IRT BCs and the Ha emission 
or the He I BCs, if all broad emission lines are sup- 
posed to share a common forming region. Johns- 
KruU & Basri (1997), analysing DF Tau, confirm 
the lack of correlation found between the IRT BC 
and the Ha emission but they found a correlation 
between the He I BC and the IRT that neither 
we nor MuzeroUe et al. (1998a) confirm in our 
samples of various CTTSs. 

We showed that the IRT lines and Ha are good 
mass accretion indicators in our sample. We also 
confirmed that the same mechanism that pow- 
ers the accretion also seems to influence the out- 
flow, as the emission and absorption components 



of Ha are well correlated with the mass accretion 
rates and with each other. Other good evidence 
for the accretion- wind connection is that the Fe 
lines in strong emission seem to be related to both 
processes, occurring when the wind and the veiling 
are prominent. 

The blueshifted absorption components are 
thought to be formed in low density winds. There 
is strong evidence that the NCs are formed at the 
stellar surface (Batalha et al. 1996). Although 
they are both influenced by the accretion process, 
they should not be compared to the theoretical 
magnetospheric profiles. 

Some of the observations tend to support the 
fact that the accretion affects the relatively quiet 
stellar atmosphere, and lines produced in that re- 
gion suffer the influence of such changes. The 
models also predict that accretion and outflow 
should be related to each other - the former occur- 
ring through closed magnetic field lines that con- 
nect the disk to the star and the latter through 
the open field lines escaping from the disk itself. 

However, the results presented in the previous 
sections also show that some of the general char- 
acteristics of the magnetospheric accretion models 
are not always conflrmed by our data. The idea 
that the broad emission components mainly arise 
from a common emitting region, the infalling mag- 
netospheric gas according to those models (Hart- 
mann et al. 1994), must be taken with care, as 
the observational evidence is not yet conclusive. 
We must also keep in mind that the magneto- 
spheric predictions do not do particularly well for 
the broad emission components that they primar- 
ily apply to. We find observed profiles to be gen- 
erally fairly symmetric with centroids distributed 
through a range of blueshifted to redshifted veloc- 
ities. Sometimes they have redshifted absorption 
components, but often not at free-fall velocities. 
Only about 20% of the BC profiles support the 
specific predictions of the models. 

Hartmann et al. (1994), MuzeroUe et al. 
(1998b) and MuzeroUe et al. (1998a) presented 
line profiles generated with magnetospheric mod- 
els and compared them with observational ones. 
They obtained good agreement between the ob- 
served and theoretical profiles of many lines for 
some stars, like BP Tau. But comparing their BP 
Tau observed profiles to ours, many differences can 
be seen. This is expected, as CTTSs are known 
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to exhibit line profile variability sometimes even 
within a few hours (Johns & Basri 1995a), and 
Gullbring et al. (1996) showed that this is the 
case for BP Tau, a typical CTTS. The problem 
that arises with the rapid profile variation is that, 
even if the magnetospheric models can reproduce 
one of the profile types of a star, it does not mean 
that all the profiles that a line may exhibit will be 
reproduced without changing basic model charac- 
teristics. As an example, the model used to com- 
pare the observed and theoretical profiles of BP 
Tau cited above does not fit our observations of 
that star well. 

Some years ago, Shu et al. (1994) proposed a 
new version of their X-celerator mechanism whose 
natural period was the rotational period of the 
star. Lines formed in magnetospheric accretion 
flows that are controlled by the stellar field might 
also show periodicity at the stellar rotation period 
in the presence of an inclined dipole field. In this 
case, the line variations would certainly be due 
to geometrical effects and a model with the same 
basic parameters should be able to describe the 
various line profiles. SU Aur presented this kind 
of variations in H/3 and Ha (Johns & Basri 1995b), 
but unfortunately its profiles do not resemble the 
published theoretical ones very much. 

In general, however, the line variations do not 
seem to be related to orbital motion as they do not 
usually correlate with the stellar rotational period. 
Different parts of a line may vary differently as 
Johns & Basri (1995a) showed by analysing time 
variations of Ha for several CTTSs. They suggest 
the line emission region is composed of discrete, 
stochastically varying blobs with a range of veloc- 
ities and turbulence, which the Sobolev treatment 
used in the actual magnetospheric models does not 
take into account. 

Najita et al. (1996), analysing Br7 emission 
profiles, showed that WL 16, an obscured low- 
luminosity YSO, presented profiles that were re- 
markably similar to the theoretical Balmer line 
calculations by Hartmann et al. (1994). Later, 
Muzerolle et al. (1998b) calculated profiles for the 
Br7 emission and showed that, in fact, the WL 16 
Br7 emission line was very well reproduced. How- 
ever, some of the stars in the sample of Najita et 
al. (1996) also present very symmetric emission 
profiles, like AS 353, DG Tau and SVS 13, that 
cannot be explained only by magnetospheric ac- 



cretion. Although the magnetospheric model can 
generically explain part of the large line widths 
and the occasional redshifted line absorptions, it 
does not predict symmetric profiles becaiise the 
magnetospheric infalling material is subjected to 
occultation effects by the star-disk system. 

Most of our broad emission profiles were found 
to be symmetric and the few asymmetries occur 
on both the blue and the red side. Due to these 
differences, the profiles presented in the theoret- 
ical papers generally do not match our observed 
ones. Many authors have speculated on the ori- 
gin of the rather symmetric broad line wings in 
TTS profiles. Basri (1990) suggested that the 
wings were formed in a turbulent region near the 
star and Edwards et al. (1994) pointed out that 
Alfven waves could generate the necessary turbu- 
lent broadening. Johns & Basri (1995b) managed 
to reproduce some of the Balmer line symmetric 
features of SU Aur by adding a range of turbulent 
high velocity components at the base of spheri- 
cally symmetric wind models and Johns & Basri 
(1995a) could also fit the DF Tau wings with a 
similar procedure. 

Most of the stars discussed here rotate slowly 
(w ^ 10km s""'^) at velocities that are much smaller 
than the velocities of the infalling material {v ~ 
200 — 300km s~^) and the inclusion of rotation in 
the theoretical models will not significantly alter 
the profiles in this case. Even in the case of a 
stiff magnetosphere, material coming from a disk 
at 3i?* would have v ~ 30km s^^ and would not 
influence the entire line. The addition of winds 
to the theoretical models will certainly change the 
line shapes but stars like AS 209 and DE Tau, 
that do not always exhibit blueshifted absorption 
and are rather slow rotators, will still need some 
further explanation of their broad symmetric line 
profiles (Figure 1). 

Stark broadening may explain the low very 
broad far wings in Ha (Hartmann et al. 1994; 
Muzerolle et al. 1998b), but will not produce the 
main Gaussian shape commonly present in many 
different lines. A better treatment of the radiative 
transfer (actually done with the Sobolev method) 
could also alter the theoretical line profiles and 
would properly take into account Stark and opac- 
ity broadening effects. 

Magnetospheric accretion probably does occur 
in CTTS. Evidence for it includes hot spots at the 
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stellar surface of these stars, thought to be due to 
the magnetic accretion shock, redshifted absorp- 
tion components at free-fall velocities that would 
be hard to explain without magnetospheric infall, 
as well as cases like SU Aur and WL 16 that give 
more explicit support to the theory. However, it 
may not be the only important mechanism. The 
symmetry of the line profiles and the lack of corre- 
lation of the variation among different line regions 
suggest an important turbulent process may be oc- 
curring together with the magnetospheric infall. 

6. Conclusions 

We have presented the spectral analysis of a 
sample of CTTSs covering a wide range of opti- 
cal wavelengths and veilings. We tested the pre- 
dictions of the magnetospheric accretion model 
and previously published results supporting them. 
We confirm that many CTTSs exhibit central or 
blueshifted broad emission lines but a substantial 
minority show redshifted centroids. The accre- 
tion and outflow seem to be related and redshifted 
absorption components are sometimes present in 
the spectra. We found that the analysis and in- 
terpretation of these is not straightforward. We 
could not reproduce the observational results that 
showed a very high frequency of redshifted absorp- 
tion in the lines of CTTSs. We also show that 
most of our emission components are symmetric, 
instead of blueward asymmetric as predicted by 
the theory or suggested by previous observational 
studies. We saw that the broad components of dif- 
ferent atomic lines are not strongly correlated to 
each other, while we would expect them to be if 
they were all formed in the accretion flow. 

Our results do not refute the general magneto- 
spheric accretion scenario, but rather indicate that 
is only part of the important processes which pro- 
duce the strong permitted line emission in CTTSs. 
Our intent is to point out that the case has not yet 
been strongly made, and that much more obser- 
vational and theoretical work is required. Part of 
the emission line profiles may also be produced in 
winds, and rotation and turbulence must be added 
to the models. 

Synoptic observations are critical to under- 
standing the full range of profiles produced by 
a single star, where some underlying parameters 
do not change but others do, and many structures 



are dynamic. 

This research is based on data collected on the 
Shane 3m telescope at Lick Observatory run by 
the University of California. We would like to 
thank Christopher Johns-KruU, Anthony Misch, 
Michael Bisset, Claude Bcrtout, Natalie Stout- 
Batalha and Celso Batalha who helped gather the 
profiles presented here. S.H.P.A. acknowledges 
support from the Conselho Nacional de Desen- 
volvimento Cientifico e Tecnologico - Brazil. 



14 



REFERENCES 

Appcnzellcr, I., Rcitermann, A., Stahl, O. 1988, 
PASP, 100, 815 

Basri, G. & Bertout, C. 1989, ApJ, 341, 340 

Basri, G. 1990, Mem. Soc. Astron. Ital., 61, 707 

Basri, G. & Batalha, C.C. 1990, ApJ, 363, 654 

Basri, G., Johns-Krull, C.M., Mathieu, R.D. 1997, 

AJ, 114, 781 

Batalha, C.C, Stout-Batalha, N.M., Basri, G., 
Terra, M.A.O. 1996, Ap.JS, 103, 211 

Beristain, G., Edwards, S., Kwan, J. 1998, ApJ, 
499, 828 

Bertout, C., Basri, G., Bouvier, J. 1988, ApJ, 330, 
350 

Calvet, N., Basri, G., Kuhi, L.V. 1984, ApJ, 277, 
725 

Calvet, N., Basri, G., Imhoff, C.L., Giampapa, 
M.S. 1985, ApJ, 293, 575 

Choi, P.I. & Herbst, W. 1996, AJ, 111, 283 

Edwards, S., Strom, S.E., Hartigan, P., Strom, 
K.M., Hillenbrand, L.A., Herbst, W., Attridge, 
J., McrriU, K.M., Probst, R., Gatley, I. 1993, 

AJ, 106, 372 

Edwards, S., Hartigan, P., Ghandour, L., An- 
drulis, C. 1994, AJ, 108, 1056 



GuUbring, E., Petrov, P.P., Ilyin, I., Tuominen, I., 
Gahm, G.F., Loden, K. 1996, A&A, 314, 835 

Hamann, F. & Persson, S.E. 1992, ApJS, 82, 247 

Hartigan, P., Edwards, S., Ghandour, L. 1995, 
ApJ, 452, 736 

Hartmann, L., Edwards, S., Avrett, E. 1982, ApJ, 
261, 279 

Hartmann, L., Calvet, N., Avrett, E., Loesler, R. 

1990, ApJ, 349, 168 

Hartmann, L., Hewett, R., Calvet, N., 1994, ApJ, 
426, 669 

Hartmann, L., Calvet, N., GuUbring, E., 
D'Alessio, P. 1998, ApJ, 495, 385 



Hartmann, L., 1998, Accretion Processes in Star 
Formation (Cambridge University Press) 

Herbig, G.H. & Bell, K.R. 1988, Lick Observatory 
Bulletin 1111. 

Herbst, W., Rhode, K., Hillenbrand, L.A. 1999, 
AAS 194, 68.03 

Johns, CM. & Basri, G. 1995a, AJ, 109, 2800 

Johns, CM. & Basri, G. 1995b, ApJ, 449, 341 

Johns-Krull, CM. & Basri, G. 1997, ApJ, 474, 
433 

Mundt, R., Eisloffel, J. 1998, AJ, 116, 860 

MuzeroUe J., Hartmann L., Calvet N. 1998a, AJ, 
116, 455 

MuzeroUe J., Calvet N., Hartmann L. 1998b, ApJ, 
492, 743 

Najita, J., Carr, J., Tokunaga, A.T. 1996, ApJ, 

456, 292 

Natta, A., Giovanardi, C. 1990, ApJ, 356, 646 

Shu, F., Najita, J., Ostriker, E., Wilkin, F., 
Ruden, S., Lizano, S. 1994, ApJ, 429, 781 

Stassun, K.G., Mathieu, R.D., Mazeh, T., Vrba, 
F.J. 1999, AJ, 117, 2941 

Valenti, J.A. 1994, Ph.D. Thesis, Univ. of Califor- 
nia, Berkeley 

Vogt, S.S. 1987, PASP, 99, 1214 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.0. 



15 



16 




Fig. 1. — Line profiles - all the data have been continuum normalized, binned and smoothed with a window 
width of 3. The horizontal dotted lines represent the continuum level at 1. The vertical dotted lines are the 
spectral line's center at the stellar rest frame. The contamination of the Na D lines by the vapor lamps of 
the city of San Jose has been removed. The dashed lines close to the IRT profiles are the velocity position 
of the blended Paschen emission lines. 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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Fig. 1. — Continued 
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